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We propose a method to calculate the QCD level density directly from the thermodynamic quan- 
tities obtained by lattice QCD simulations with the use of the maximum entropy method (MEM). 
Understanding QCD thermodynamics from QCD spectral properties has its own importance. Also 
it has a close connection to phenomenological analyses of the lattice data as well as experimental 
data on the basis of hadronic resonances. Our feasibility study shows that the MEM can provide 
a useful tool to study QCD level density. 
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1. QCD level density 

In this report, we discuss the QCD level density in a finite box with a size V. Denoting energy 
eigenstates of the QCD Hamiltonian J^qcd by E n , the energy-level density A(E,V) is defined by 
counting the number of states in the range of E to E + dE, i.e., 

A(£,V) = £S(£-£ n ) with Jf QCD \*¥)=E n \y). (1.1) 

n 

Using this level density, the partition function is written as its Laplace transform; 

-i f'OO 

3?(p,V)=Tr g-^ecfl = / A{E,V)e^ E dE, (1.2) 
L J Jo 

where j3 = 1 JT with T being the temperature of the system. Thermodynamic quantities such as 
the pressure, p, and the energy density, e, are related to the partition function as 

KW-^ow, «W = -i^W (13, 



Eq.(1.2) together with Eq.(1.3) imply that every thermodynamic quantity is obtained from the in- 
formation of the level density defined at zero temperature. Especially, Eq.([L2|) tells us that the 
properties of QCD phase transition at finite T, which is characterized by the non-analytic behavior 
of p(P) in the thermodynamic limit, are already encoded in the QCD level density defined at T = 0. 

The investigation of the level density in nuclear physics and in hadron physics has a long 
history: For example, Bethe studied the nuclear level density by evaluating the partition function 
of the Fermi gas [p. By performing the inverse Laplace transform, he showed that A oc exp(2\/c£ r ) 
with c being a constant related to single-particle level density. Hagedorn studied hadronic level 
density and has derived an asymptotic formula for the state density of hadrons with a mass m, 

p(m) — exp(m/r H ). (1.4) 



Here a and 7# are some constants and the latter is called the Hagedorn temperature fl2fl. Eq.(|1.4|) is 
derived from the celebrated bootstrap model in which the state density p (m) and the energy-level 
density A(E = m,Vo) in a fireball of size Vq are identified for large m. The exponential growth of 



the hadronic state density, Eq.(1.4), agrees with experimental data up to 2 GeV at present and the 
agreement becomes better as new resonances are included [pp. Because of this success, the hadron 
resonance gas model as well as the Hagedorn's formula have been and is being widely used in QCD 
phenomenology. However, this model describes only the hadronic matter at temperature below 7#. 



If we try to calculate 3f with Eq.(|L4|) above 7#, the integral does not converge. 

On the other hand, progresses in QCD thermodynamics has been obtained steadily by the first 
principle lattice simulations. The results show that e and p increase rapidly near the transition 
temperature T c and approach to the black-body formula, e ~ 3p oc T 4 , for sufficiently high temper- 
ature. It has been also found that thermodynamics below T c is well described by hadron resonance 
gas model as long as appropriate hadron masses relevant to lattice simulations are employed [Q,|5|]. 
Therefore, the time is now ripe to consider a unified and model-independent description of QCD 
thermodynamics on the basis of the QCD level density. In fact, detailed studies of the level density 
is particularly useful to identify the relevant degrees of freedom in hot QCD below and above T c . 
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Figure 1: A schematic model for the pressure (left panel) and the energy density (right panel). 



Now let us first pay attention to a close similarity between Eq.(1.2) and the spectral repre 



sentation of the hadronic two-point correlation function, G, as a function of the imaginary time 



T: 



G(z) 



R(co)e-™dCQ, 



(1.5) 



where R(oo) is the spectral function and oo is the frequency. Formal correspondences, i2° <-> G, 
j3 <-> T, E <-> oo and A <-> 7?, are clear by comparing Eq.(1.2) with Eq.(1.5). Since the maximum 



entropy method (MEM) is known to be a powerful tool to extract the spectral function R from the 
lattice data G the same technique is expected to be used to extract the level density A(E,V) 
from the lattice data of 



2. A toy model 

Before testing the idea of using MEM to extract the level density A(E,V), let us first discuss 
the general structure of A expected in simple cases where some analytic study is possible. First of 
all, the inverse Laplace transform (the Bromwich integral) of Eq.(L2) reads 

' y—i° 

where / is a real number chosen so that all the singularities of £*f(P,V) are to the left of it. In the 
leading order of the saddle point approximation, one readily finds 



1 / i 7+i°° i / i 7+/°° 

A(E,V) = — 3?{p,V)e Ep dl3 = — e MP)v+E)p dp, (2.1) 

Altl J J— /oo L%1 J y /co 



A(E,V)«exp(s(e)V), (2.2) 

where e(=E/V) is the energy density and s(e) is the entropy density. 

Consider a simplest case of free massless particles where the pressure is given by p(j3) = 
(a/3)r 4 . Then, by working out the Gaussian integral around the saddle point, one finds 

A(E,V )lEVU>>>1 ~ -±= (^) " 8 exp feWW) , (23) 
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Figure 2: Left panel: A schematic picture of s(e) lnA(.E,V) in a toy model as a function of e = E/V. 
Right panel: The QCD level density A (E, V) in an arbitrary unit extracted from the lattice data. 



which implies InA ~ Z? 3 / 4 at high excitation energies. 

To examine a more realistic example with a phase transition, let us consider the bag equation 
of state in which we assume free massless pions (free massless quarks and gluons) below (above) 
T c . In this case, the pressure of the system is given by 



P(T) 



V 



injr 



T 4 6(T c -T) + ^T 4 -Bjd(T-T c ), (2.4) 

where B is the bag constant which is related the critical temperature as T c = [(3fi) / (o q — a/j)] 1 / 4 . 
a/, and a q are proportional to the number of degrees of freedom in the hadronic phase and the 
quark-gluon plasma phase, respectively. A schematic sketch for p/T and e/J 4 is given in Fig. 1. 



For the bag equation of state given in Eq.([2i4), the exponent s(e) in Eq.(|2.2[) turns out to be 
4 



-a 1/4 £ 3 / 4 



[I]- 



1 ^3 £ 



II 



4 



3 a « 



V4, 



B 



,3/4 



HI 



(2.5) 



Here the formulas [I], [II] and [III] are valid for e < £ c \, E c \ < £ < £ C 2, and e > £ c 2, respectively. 
Note that e c \ = Of,T 4 {e c2 = o q T 4 + B = | (4o q — Of,) T 4 ) is the energy density just below (above) 
the phase transition point as shown in the right panel of Fig. 1 . 



The behavior of the exponent s(e) given in Eq^^) is illustrated in the left panel of Fig. 2. 
The level density shows the Hagedorn type behavior In A oc E in the phase transition region [II], 
while it shows softer behavior In A ~ £ 3 / 4 in the high temperature region [III]. The exponent has a 
crossover from E to Z? 3 / 4 at E = E c i = £aV at which the quark-gluon plasma starts to emerge. 



3. Application of MEM 

Let us now try to extract A(E,V) from lattice QCD data of the thermodynamic quantities. 
Following the idea of MEM, we minimize the "free-energy" functional <2(A) with respect to the 
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level density A: 



<2=y-«S inf . (3.1) 



Here, Sm is the information entropy given by 

S int = J\f-l-flnf)dE, A(E,V)=A asy (E,V)xf(E,V), (3.2) 

where A asy (E, V) implies a "default model" representing the asymptotic behavior of A in the large 
E limit. In the present study, we adopt A !iS y(E,V) = (V 1 ' 3 / VSn)exp (|a'/ 4 Ve 3 / 4 ) which is pro- 
portional to the asymptotic form of A for a free quark-gluon gas. We adjust a to reproduce £ /T 4 
measured on the lattice at high temperature. For the likelihood function # 2 , we have chosen 




^\aJA(E,V)e-^ E dE 



J{E/V)A(E,V)e-l i i E dE 



,Jat \ 



where pf x and £- at are the data obtained by lattice simulations at n discrete values of the inverse 
temperature /3,-. Also, Ap| at (Ae ; ! at ) denotes statistical error of the pressure (energy density). 

We use the 2-flavor full QCD data with improved Wilson quarks generated on a 16 3 x 4 lattice 
[R]. The spatial size of the lattice in the physical unit is approximately (4fm) 3 , although the physical 
volume changes for different values of /$;. As a first step, we choose V = 1 fm 3 in Eq.(3.3) by 
assuming small volume-dependence of the lattice data. We use the data for the pressure and energy 
density with mps/my = 0.90 in which there are "seven" independent data points (n = 7). 

A(E,V) reconstructed by MEM is shown by the solid line in the right panel of Fig. 2 in an 
arbitrary unit with the logarithmic scale. Dashed line is A asy mentioned above. To set the scale in 
the horizontal axis, we assume T c « 175MeV. Also, we have chosen a = 1 in the present MEM 
analysis: eventually it has to be eliminated by calculating the probability distribution P[a] [gfl. As 
can be seen from the figure, the behavior of A(E,V) at high E is consistent with A 

asy rs "' 

£ 3 / 4 (the 

dashed line), while A(E,V) decreases strongly around E = 9 — lOGeV and deviates substantially 
from A asy at low energies. This rapid crossover of InA is qualitatively consistent with what we have 
discussed using the toy model. 

Now we briefly discuss some systematic uncertainty in the MEM analysis. Instead of utilizing 



both pf l and £- at as in Eq.flOp, we have done MEM analyses by using pf l only and by using e' at 



only. In these cases, the crossover region has moved ±20% from that shown in the right panel 
of Fig. 2. Therefore, the systematic error due to different choices of % 2 is still large at present as 
long as we take only seven data points. Nevertheless, we believe that MEM could become a useful 
tool to extract the QCD level density if large number of data points with high accuracy become 
available in the future. 

4. Conclusions and outlook 

To understand thermodynamic properties near the QCD phase transition, the knowledge of 
the energy-level density is quite useful. In this report, we proposed a new method to calculate the 
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QCD level density and have made a feasibility test using lattice data. Although further systematic 
studies are necessary by increasing the number of data points, the direct calculation of the QCD 
level density seems to be possible using the maximum entropy method. 

An extension of the present study to the system at non-zero baryon density is interesting to 
be explores in connection with recent progress of finite density lattice QCD. The grand partition 
function 3f(f},[i,V) can be separated into the canonical partition functions 3f^(P,V) for each 
fixed quark number N by the fugacity expansion as 

iT(j3,M,V)=£^^(j3,V)=£e A '^ A(E,N,V)e~^ E dE. (4.1) 

Then, by using the lattice data with several different values of j3 and jJ,, one may extract the level 
density A(E,N,V). Such an analysis will shed lights on colored composite states with non-zero 
baryon numbers above T c such as the quark-gluon bound state and the diquark [|| ||]. If such states 
are important to thermodynamic quantities, they should also show up in the QCD level densities 
with N = 1 and N = 2. 
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